Dynamic and spectral mixing in nanosystems. 
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In the framework of simple spin-boson Hamiltonian we study an interplay between dynamic and 
spectral roots to stochastic-like behavior. The Hamiltonian describes an initial vibrational state 
coupled to discrete dense spectrum reservoir. The reservoir states are formed by three sequences 
with rationally independent periodicities 1 ; 1 ± 8 typical for vibrational states in many nanosize 
systems (e.g., large molecules containing CH2 fragment chains, or carbon nanotubes). We show 
that quantum evolution of the system is determined by a dimensionless parameter ST, where F 
is characteristic number of the reservoir states relevant for the initial vibrational level dynamics. 
When 5 T > 1 spectral chaos destroys recurrence cycles and the system state evolution is stochastic- 
like. In the opposite limit 8T < 1 dynamics is regular up to the critical recurrence cycle k c and 
for larger k > fc c dynamic mixing leads to quasi-stochastic time evolution. Our semi-quantitative 
analytic results are confirmed by numerical solution of the equation of motion. We anticipate that 
both kinds of stochastic-like behavior (namely, due to spectral mixing and recurrence cycle dynamic 
mixing) can be observed by femtosecond spectroscopy methods in nanosystems in the spectral 
window 10 11 - 10 13 s" 1 . 
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Quantum dynamics for an initial vibrational level e?coupled to a set of discrete dense levels e° can be described in 
a framework of so-called spin-boson Hamiltonian [l[ , [2[ 

H = e° s b+b s + Y,£nKK + + b nK) > (1) 

n n 

where 6+(6 s ) is initial level creation (annihilation) operator (i.e., excitation of the initial vibrational level from the 
system ground state, which is assumed so deep that its influence on the system dynamics can be neglected), and 
b^(b n ) are similar operators for the discrete reservoir levels. C' n stands for the coupling matrix elements. For this 
Hamiltonian time dependent wave function can be presented in a series over full orthogonal set of wave functions of 
the unperturbed (uncoupled) initial and reservoir states with time dependent coefficients (amplitudes) a s (t) 




a s(t) = l^il^) eM-iEt)} , (2) 

E=e n 

where F(E) = is a secular equation to find the eigenstates e„ of the Hamiltonian ([TJ. 

Common wisdom Q , Q , Q , [(| claims that stochastic-like dynamics is a feature typical for random matrix Hamil- 
tonians, with random eigenstate spectra. However for a system with discrete dense spectrum (e.g., vibrational states 
in nano-particles or medium size molecules) there is another dynamic root to quasi-stochastic behavior. Indeed for 
such a system with discrete spectrum, dynamic evolution is represented by periodically repeating steps (recurrence 
cycles). When time is going on, the initial vibrational level population oscillates faster and faster and corresponding 
response signals become broader. Eventually at a certain critical cycle number k c , the cycles are overlapped in time. 
Then for any finite accuracy of time or frequency measurements, time evolution looks as irregular, quasi-stochastic, 
indistinguishable from truly chaotic behavior. It was demonstrated recently 0], @], @ for the simplest version of the 
spin-boson Hamiltonian ((J) with = n, and C n = C (so-called Zwanzig approximation [To|). Here (and in what 
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follows) we utilize a characteristic reservoir spacing as energy unit, and n = is the reservoir level coinciding with e s . 
It is worth noting that this dynamic mixing occurs in a single Hamiltonian system (not for an ensemble of equivalent 
systems). 

Since there are mentioned above two roots to stochastic-like time evolution, two natural questions arise, namely, 
(i) how the both mechanisms are interrelated, and (ii) whether they are independent ones. To answer these questions 
is a purpose of our presentation. We study quantum dynamics for a version of spin-boson Hamiltonian (|T|), where 
the both ingredients yielding to stochastic time evolution (spectral and dynamic) may be presented simultaneously 
and can be tuned by model Hamiltonian parameters. We assume that the reservoir discrete bare spectrum in the 
Hamiltonian (TTJ) is formed by three equidistant sequences with different periods 



±3n : 



±(1 - 8)(3n + 1) ; e° 2 (n) = ±(1 + 8)(3n + 2) . 



(3) 



where spectral shift dimensionless parameter 8 < 1/2. It is worth noting that such kind of triplet structures are quite 
frequently observed in the spectra of many molecular systems, forming notorious Snyder sequences 11 1| (see also more 
recent publications on high resolution spectroscopy of nanotubes [l2| , phospholipid molecules [l3[ , [15| and fullerenes 
[HJ|). When the total number of levels N — > oo and the sequence periods 1, 1± 8 are rationally independent, the 



spectrum of the reservoir becomes mixed. In own turn, according to the ergodic properties jl6j, 17| chaotic behavior 
of an ensemble of systems (random eigen- value distribution) holds if its individual system is mixing. That is why we 
term our case as exhibiting of stochastic-like behavior (cf. with definition [l8j). 



In more practical terms mixing phenomenon can be related to reservoir level ordering. Indeed, level ordering and 
interlevel spacings depend on a cycle number. In a zero cycle the spectrum is formed by triplets in the following 
sequence 



e 2 {n - 1) < e (n) < ei(n) < e 2 {n) . 



(4) 



When n increases the quantities p^ = eo{ n ) ~~ — 1) an d p% = ei(n) — eo(n) become smaller and fill with a 
step 38 intervals [0 , 1 ± 8]. Level splitting between neighboring triplets p^ = ei(n) — eo(n) increases with a step 68 
approaching to 3. In the next cycle k = 1 the level ordering is different 



ei(n) < eo(n) < e 2 (n - 1) < ei{n + 1) . 



(5) 



and the quantities Pi (n) = 62(71 — 1) — eo(^) and ?4 ( n ) — e o( n ) ~ e i{ n ) are increased with the same as in the k = 
cycle, step 38 up to a limit value 3/2. One can check that in the next even cycles the splittings Pi and p% decrease, 

whereas p^^ increases. In the odd cycles the opposite variation of the splittings hold (jpi and p 2 increase, and P3 
decreases). Eventually levels from different triplets fill in the k — > oo limit almost uniformly and densely the interval 
[0 , 68] with odd-even alternations with the recurrence cycle number. By rather lengthy and boring but straitforward 
calculations we find that for cycle number k > 1, the variables p\ and p 2 fill uniformly the intervals 
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with the different order of p\ and p 2 in even and odd cycles. The P3 fills the interval 

3 
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(7) 



In own turn the above intervals limits a 1 with j = 1,2,3 satisfy the following relations 

.<«_u (3t+2)± a ; „<»» = a ( 3 t+2) _i 



(8) 



where {X} stands for the fractional part of X. At k ^> 1 the (jSJ can be generated approximately by so-called 
fractional part recurrence relation 



M 

l k+l 



28 Qk 



(9) 



As it is known [19[ the fractional part transformation is the mixing one (i.e., stochastic) with characteristic correlation 
time (in our notations) oc ln(l/<5). Thus our model reservoir spectrum is the reservoir with spectral mixing. This 
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phenomenon can be formulated as the mixing of pj parameters within the interval [0 , 3] . The interlevel spacing 
distribution density is approximately (over the parameter 1/iV) a constant ~ 5/9 in the interval [0 , 3/2] and another 
constant ~ 1/9 in the interval [3/2, 3]. Without going to more subtle mathematical details of random sequences 
we calculate numerically the distribution function p(e) for our model spectrum. When the total number N of the 
reservoir levels increases the distribution approaches to that with two uniformly distributed parts (see the fig. 1, 
where we show also the widely used in the literature Wigner distribution (see e.g., 0,01, and also [20] ) which holds 
for orthogonal Gaussian random matrices |5j). 
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Figure 1: Interlevel spacing distribution density p(e) for the reservoir ([3]) with 3 rationally independent periods (S — 0.049). 
Empty circles (C 2 = 0), stars (C 2 = 1), and filled circles (both distribution densities coincide) are density 5/9 in the interval 
[0 , 3/2] and another constant 1/9 in the interval [3/2 , 3]. Wigner distribution is shown for comparison by dashed line. 

Let us come back to the secular equation to find the eigenstates 

n=+N , „ 2 C' 2 C 2 \ 

F{e) = e- £ (^ v + ^+i_ + ^+^)=0. (10) 

n= _ w V e £ 3n 6 e 3ri+l 6 £ 3n+2 / 

For the ease of algebra (if needed this approximation can be relaxed) we assume following Zwanzig 10] approximation 

Csn = C3„+i = Czn+2 = C 

and in this case the series entering (|10j) can be written explicitly in terms of trigonometric functions 

Here T = irC 2 characterizes the window of the reservoir states, where at S = the coupling to the initial vibrational 
level e s = contributes essentially to system time evolution (the reservoir levels with the quantum numbers < T 
govern the system dynamics) . 

When S — our model is reduced to a single equidistant reservoir Zwanzig model, and in this case dynamic 
stochastic-like behavior (due to recurrence cycle mixing) occurs 0, @ for the cycle number k > k° c = ttT. If S + 
analysis is more involved. One can check by tedious but direct calculations that the trigonometric equation (|11[) has 
one and only one root in each interval between the values of the bare energies e^n) , t\{n), and e% (n) given by the 
expressions ([3]). Armed with this knowledge the formal solution @ for the initial vibrational level amplitude a s (t) 
can be represented as a finite series over the partial recurrence cycle amplitudes fits (r/s) (where k is confined within 
the interval < k < [t/(2ir)] and [y] means an integer part of y) calculated in the local cycle time = t — 2irk. 
In own turn ai (rk) can be presented in terms of generalized Lommel and Laguerre polynomials |2l| . To avoid this 
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lengthy and cumbersome mathematics in this publication we restrict ourselves to the only small 5 <C 1 limit, when 
the secular equation is reduced to 



e-rcot(7re)[l + R(e , 5)} = . 



(12) 



with 



R(e,5)= [l + q 



sin(7re/3) 
sin(7re) 



-l r 



sin(7re/3) cos(7re/3) 
sin(7re) 3cos(7re) 



„sin(7re/3) , r- 
-25 V '> (V3 + q) 
cos(7rej 



(13) 



where 



<7 = 1 + 4sin(7ri5e) cos(7rJe + 7r/3) . (14) 
All qualitative features of the quantum time evolution for our model can be seen in this limit (5^1. Since F(e) 

(k) 

depends on ire/ 3 in each recurrence cycle k its amplitude a s has not only the main Loschmidt echo signal but also 
two sattelites shifted with respect to the main signal by 27r/3 and 47r/3 (see fig. 2 where we plot initial vibrational 
level population time evolution). Intensities of the sattelites are proportional to 5, and when the cycle number k 
increases the main echo signal and its sattelites acquire a fine internal structure. The number of components in the 
fine structure of the echo signals increases with k and starting with a certain threshold k c start to overlap. For 6 <C 1 

fc" 1 = (tiT)- 1 +3<5. (15) 

At k > k c the triplets formed by the main echo signal and sattelite triplets are strongly mixed and time evolution 
becomes chaotic. In the fig. 3 we computed k c (S). It is clear that the numerical results are in a good agreement 
with our approximated formula (|15|) . This formula (|15l) is our main result in this paper. It shows how the spectral 
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Figure 2: Initial vibrational level population dynamics (C 2 = 1). From the top to the bottom: 5 — ; 8 = 0.019 ; S = 0.049 ; <5 = 
0.079 ; S — 0.12. Critical cycle numbers k c are indicated by arrows. 

source of the chaotic behavior (6 in the (|15p ) interplays with the dynamic source of stochastic- like time evolution (r 
in the (|15I0 . Namely, k c {5) is determined by the condition that the triplets from the cycle k c are mixed with those 
from the next cycle. The spectral chaos contribution (5 ^ 0) reduces considerably the number of cycles with regular 
dynamics. For S = 0.12 (still 5 <C 1) the intensities of sattelites and the main echo component are comparable and 
the regular triplet structure is completely broken after 2 initial cycles. We illustrate this in the fig. 2, where we show 
numerically calculated |a s (£)| 2 (i.e., the initial vibrational level population) for 5 — , 0.019, 0.049, 0.079, 0.12 (a 
bit bizarre numerical values of S are chosen to get better approximation for rationally independent sequences of the 
reservoir levels). 

Thus we arrive at the general conclusion. Dynamic, stochastic-like time evolution is determined by the local 
spectrum characteristics (mean interlevel spacing) and occurs only for sufficiently large recurrence cycle number. On 
the contrary spectral truly chaotic dynamics is governed by the global spectrum structure (e.g., in our spin-boson 
model Hamiltonian (JTJ) by the reservoir triplet structure ([3])). It is instructive also to show contributions of the 
reservoir states into the initial vibrational level wave function amplitude a s (t). According to the expression @ we 
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Figure 3: Critical cycle number dependence on 8: C 2 = 1 - open circles; C 2 = 2 - stars; C 2 = 4 - squares. 



computed the quantity A n = (dF/dE)~ 1 \E=e n and plotted the results in the fig. 4. These quantities A n also manifest 
themselves simultaneous effects of dynamic and spectral mixing. For the pure Zwanzig model (no spectral mixing, 
i.e., (5 = 0), the coefficients A n decrease monotonously with e„. However, upon increasing 5 (spectral mixing), A n 
distribution becomes more and more irregular. 
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Figure 4: Amplitudes A n = (dF/dE)' 1 ^ © of the reservoir states contributing into the initial vibrational level evolution. 
From the top to the bottom: 5 = ; 5 = 0.049 ; 5 = 0.12, and C 2 = f . 

Our motivation in this work is not a pure curiosity. As a matter of fact quantum dynamics of various systems 
(ranging from relatively small molecules up to large photochromic molecules and their protein complexes or molecules 
confined near interfaces [22| (see also 23]) is an active area of experimental researches. The femtosecond spectroscopy 
data (which allow to study time evolution of one initially prepared by optical pumping state) manifest variety of 
possible regimes including not only weakly damped more or less regular oscillations but also very irregular long 
time behavior with a number of peaks corresponding to a partial recovering of the initial state population. This 
generic feature is omnipresent in the systems with complex and irregular vibrational relaxation. The triplet model 
investigated in this paper reflects the spirit of minimalist approaches, in that it is simple yet based on a physical 
principle. Understanding all its limitations, we nevertheless hope that our crude theory captures the essential elements 



6 

of vibrational relaxation in nano-systems. 
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Figure captions 

Fig. 1 

Interlevel spacing distribution density p(e) for the reservoir ^ with 3 rationally independent periods (S — 0.049). 
Empty circles (C 2 = 0), stars (C 2 = 1), and filled circles (both distribution densities coincide) are numerical compu- 
tation performed for 166 levels (10 cycles). Solid line - our theoretical prediction: p(e) ~ 5/9 in the interval [0, 3/2] 
and p(e) ~ 1/9 in the interval [3/2 , 3]. Wigner distribution is shown for comparison by dashed line. 

Fig. 2 

Initial vibrational level population dynamics (C 2 = 1). From the top to the bottom: 6 = 0; 6 = 0.019; 5 = 
0.049 ; 5 — 0.079 ; 5 = 0.12. Critical cycle numbers k c are indicated by arrows. 

Fig. 3 

Critical cycle number dependence on i5: C 2 = 1 - open circles; C 2 = 2 - stars; C 2 = 4 - squares. 
Fig. 4 

Amplitudes A n = (dF/dE)~L e © of the reservoir states contributing into the initial vibrational level evolution. 
From the top to the bottom: 8 = ; S = 0.049 ; S = 0.12, and C 2 = 1. 



